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In this work, the charge transfer fluctuation which was previously used for pp collisions is proposed 
for relativistic heavy-ion collisions as a QGP probe. We propose the appearance of a local minimum 
at midrapidity for the charge transfer fluctuation as a signal for a QGP. Within a two-component 
neutral cluster model, we demonstrate that the charge transfer fluctuation can detect the presence of 
a QGP as well as the size of the QGP in the rapidity space. We also show that the forward-backward 
correlation of multiplicity can be a similarly good measure of the presence of a QGP. Further, we 
show that the previously proposed net charge fluctuation is sensitive to the existence of the second 
phase only if the QGP phase occupies a large portion of the available rapidity space. 

PACS numbers: 24.60.-k, 25.75.-q, 12.38.Mh. 


I. INTRODUCTION 


Active researches in relativistic energy heavy-ion col¬ 
lisions have given us much information about the hot 
matter produced in such collisions. Much attention has 
been directed to the question of whether a deconfined 
quark-gluon plasma (QGP) phase has been formed. The 
experimental studies do suggest that strongly interacting 
dense matter was formed during the early stage of reac¬ 
tion, and the energy density of such matter is very high 
(see Hiiii and references therein). Most theoreti¬ 
cal models for such hot and dense matter explicitly invoke 
quark andg^on degrees of freedom in the elemen tary 
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One way to detect the presence of a QGP is then to 
measure the changes in the fluctuations and correlations 
which could originate from the new phase of matter. 

In this work, we propose charge transfer fluctuations 
as a signal of the presence of a QGP as well as the mea¬ 
sure of the (longitudinal) size of the QGP. The charge 
transfer fluctuation for elementary collisions was origi¬ 
nally proposed by Quigg and Thomas where they 
considered a flat charged particle distribution dNch/dy. 
This idea was later extended to smooth distributions by 
Chao and Quiggplj. 

The central result of Refs. min is the relationship 
between the single particle distribution function dNch/dy 
and the charge transfer fluctuation: 


Du{y) 


dNch 

dy 


where k is a constant and 
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is the charge transfer fluctuation. The charge transfer 
u{y) is defined by the forward-backward charge differ¬ 
ence: 


u{y) = [QF{y)-QB{y)]/2, (3) 

where Qp (y) is the net charge in the rapidity region for¬ 
ward of y and Quid) is the net charge in the rapidity 
region backward of y. The fluctuation Du{y) is then a 
measure of the correlation between the charges in the 
forward and the backward regions separated by y. 

The importance of the relationship o lies in the fact 
that K is in fact directly proportional to the local unlike- 
sign charge correlation length. Heuristically, this can be 
explained in the following way. Suppose all final parti¬ 
cles originate from neutral clusters and each cluster pro¬ 
duces one positively charged particle and one negatively 
charged particle. Then the only way u{y) can deviate 
from zero is when one charged particle from a cluster 
ends up in the forward region while the other ends up in 
the backward region as illustrated in FigHi For each one 
of these split pairs, the charge transfer u{y) undergoes 
a 1-D random walk with a step size 1. Therefore, the 
charge transfer fluctuation Du{y) should be proportional 
to the number of split pairs, or equivalently the number 
of random steps taken. 

If A is the typical rapidity difference between the two 
decay particles from a single cluster, then only the clus¬ 
ters within the rapidity interval {y — Xj2,y + A/2) can 
contribute to Du{y) as illustrated in FigQ] The num¬ 
ber of such clusters is then XdNdstr/ dy where dNc\str/dy 
is the density of the clusters at y. Since the final 
particle spectrum dNch/dy should be proportional to 
dNc\str/dy, we have Eq.(^ with k oc A. Hence the ratio 
K = Du{y)/{dNch/dy) is a measure of the local environ¬ 
ment near y: If A is a function of y, then k(j/) should also 
change accordingly. 

The fact that k is constant in elementary collisions 
indicates that in such collisions the correlation length 
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FIG. 1: A schematic illustration of the charge transfer fluc¬ 
tuations in the rapidity space. Only the pairs within A/2 of y 
can contribute to the charge transfer fluctuation Du{y)- Here 
A is the rapidity correlation length, or the rapidity distance 
of the decay particles from a single cluster. If A is a function 
of y, then Du{y) also changes with y. 


is constant throughout the entire rapidity range (see 
mill and references therein). However, if a QGP is 
produced in the central region of the relativistic heavy 
ion collisions, we can expect the local charge correlation 
length 7 ( 2 /) increases as y moves away from central ra¬ 
pidity. This is because the charge correlation length in 
a QGP is expected to be much smaller than that in a 
hadronic phase |3 ■ Iii this case, the ratio 




(4) 


will vary from a smaller value to a larger value as one 
goes away from the central region toward the forward 
region. 

There have been many studies of the fluctuations and 
correlations in heavy ion collisions 1 ^ 1 ^ 

[13, il 111 111 HIM El mill III ifl EOT Most 

of these studies concentrate on global information and do 
not address possible spatial inhomogeneity of the created 
matter in relativistic heavy-ion collisions. For instance, if 
the QGP phase is confined to a small rapidity region, the 
net charge fluctuation measures proposed in Refs.|13IU 
may not be very sensitive to the presence of the QGP. 
Hence negative results from experiments I11I13E3. do 
not necessarily exclude the formation of a QGP. 

Our expectation that the central rapidity region in 
the heavy ion collisions is mostly QGP originates from 
Bjorken’s seminal work [H. In that paper it was as¬ 
sumed that the expanding QGP evolves in a boost- 
invariant manner. Such an assumption naturally leads 
to the expectation that the central plateau in the ra¬ 
pidity spectrum is a manifestation of a boost-invariant 
QGP. However, recent RHIC results cast some doubts 
on boost-invariant scenario in the central rapidity region: 
Although the charged particle distributions as a func¬ 
tion of pseudo-rapidity shows a central plateau 0 ,E 3 l, 
the recent rapidity spectrum of charged particles from 
the BRAHMS group is consistent with a gaussian follow¬ 
ing the Landau picture |4 (t| although a plateau within 
— 1 < y < 1 cannot be ruled out ^3- The elliptic 
flow spectrum from the PHOBOS group ^3 shows no 


discernible plateau at all as a function of the pseudo¬ 
rapidity. Thus, a simple boost invariance scenario in 
a large ran ge o f rapidity space as originally envisioned 
by Bjorken may not be valid. If the QGP phase is 
produced in the relativistic heavy-ion collisions, a pure 
phase may very well be confined to a very limited rapid¬ 
ity range. It is, therefore, important to have an observ¬ 
able that is sensitive to the local presence of a QGP. The 
charge transfer fluctuation is such a local measure of a 
phase change. 

Of course, as emphasized in Ref. ^3li ^ particular type 
of fluctuations is just one particular aspect of the under¬ 
lying correlations. Usefulness of each type of fluctuations 
then depends on the sensitivity of the chosen fluctuation 
to an interesting aspect of the correlation. For charge 
transfer fluctuations, that aspect is the size of the local 
charge correlation length. Hence if the QGP phase is spa¬ 
tially confined to a narrow region around the midrapidity, 
the charge transfer fluctuations can signal its presence 
and also can yield information about the size. 

In this study, we propose the appearance of a clear 
minimum at midrapidity for the ratio K{y) as a signal 
for the existence of two different phases. The slope and 
the size of the dip around midrapidity can then reveal 
the size of the new phase (presumably a QGP). These 
features should disappear as the energy is lowered or the 
collisions become more peripheral where a QGP is not 
expected to form. 

In the following, we use a single component neutral 
cluster model and a two component neutral cluster model 
to study the purely hadronic case and the mixed phase 
case. However, the fact that the charge transfer fluctua¬ 
tion is a useful measure of the local correlation length is 
independent of our particular choice of models. Hence we 
expect that the general conclusions drawn in this study 
should be valid even within more sophisticated models as 
well as in real experimental situations. A case study us¬ 
ing cascade models with an embedded QGP component 
is under way. 

We note here that most of the discussions in this study 
are in terms of the rapidity y. However, the validity 
of our results does not depend very much on whether 
rapidity y is used or the pseudo-rapidity 77 is used. We 
also note here that the argument given here applies with 
very little change to any conserved charges such as the 
baryon number. 

The rest of this paper is organized as follows: In the 
next section, we consider the basic phenomenology of the 
charge transfer fluctuations. In Sect IHIl we consider the 
net charge fluctuations and the charge transfer fluctua¬ 
tions in a single component model. In Sect llVI we present 
our main results on a two component model. It is pro¬ 
posed that the presence of a rising segment of the charge 
transfer fluctuation as a function of rapidity can be used 
as a QGP signal. We also show that the charge transfer 
fluctuation can reveal the size of the QGP. A summary 
is given in Sect El 
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FIG. 2: The charge transfer fluctuation results from pp col¬ 
lisions at the beam energy of 205 GeV is shown as a function 
of rapidity. The line is the charge yield profile measured in 
the same reaction scaled by a factor of 0.62. This is a re-plot 
the results reported in 1^ . 


II. CHARGE TRANSFER FLUCTUATIONS 


The charge transfer is defined in Eq. 0 . The charge 
transfer fluctuation is defined in Eq. Originally Quigg 
and Thomas |3, considering a flat dNch/dy, argued that 
if all hadrons originated from neutral clusters, then the 
following relation should hold: 


^ max 


(5) 


where A is the rapidity correlation length of unlike-sign 
(H—) pairs originating from a single neutral cluster. Nch 
is the total multiplicity of the produced charged particles 
and is the beam rapidity in the CM frame. 

Later, this was extended by Chao and Quigg to 
smooth charged particle distributions to yield Eq. o, 
Du{y) = KdNch/dy, with k oc A. The experimental re¬ 
sults on pp and K~p collisions show that this relationship 
is remarkably good with k ranging from 0.62 to 0.85 (see 
mill and references therein). We have re-plotted 205 
GeV pp collision results from Ref.l^^ in FigH Here the 
proportionality constant k is approximately 0.62. As we 
shall see later in this section, the proportionality of the 
charge transfer fluctuation to the charged particle spec¬ 
trum is a strong argument for correlated charge pairs 
instead of uncorrelated charged particles. 

The charge correlation length A measures the rapidity 
correlations between unlike-sign charges. (This quantity 
also plays a central role in two of the proposed QGP sig¬ 
nals, namely the net charge fluctuations and the balance 
function.) To illustrate the relationship between the pro¬ 
portionality constant k and the unlike-sign charge cor¬ 
relation length A in Eq. 0, we consider a simple ‘p’ gas 
model (see e.g. |5?i|l. 

In this model, each is assumed to decay to a 
pair and each ‘p^’ is assumed to decay to a tt+tt” pair. 


This is similar to the ‘p’ and ‘w’ models used for pp colli¬ 
sions [^1^. One should not, of course, regard these p’s 
as physical p-mesons. These are just convenient names 
for charged and neutral clusters. In particular, they are 
not isospin triplets. 

Consider a set of events where Mq number of p°’s and 
M+ and M_ number of p+’s and p“’s are produced. The 
full joint probability for the rapidity of the charged pions 
for this set is given by: 

M+ M_ Mo 

piiva}) = Y[9{yi)Y[9{yj) n /o(yfc> ( 6 ) 

i=l j—1 k—1 

where fo{y^,y~) is the probability for the two decay 
products of a neutral cluster to have the rapidities y^ and 
y~ and g{y) is the single particle distribution function 
for the charged particles originating from the charged p’s. 
Averaging over the distributions of Mq, M_|_, M_ with the 
condition Q = M+ — M_ = constant, it is not hard to 
show that 

/ y /-oo 

dy' / dy” foiy',y'') 

-oo J y 

/ y poo 

dy'yiy') / dy" g[y”). (7) 

-OO J y 

where Mch = Af_|_ -|- M_. In arriving at the above re¬ 
sult, we neglected a contribution that is on the order of 
{u{y))‘^/Nch ~ Q^/Nch where Wh = 2(Mo)-b(Mch). This 
should be small compared to the terms in Eq. 0 when 
Q <C Nch- See Appendix B for details. 

The single particle rapidity distribution is given by 

= 2 (Mo) h{y) -b (Mch) p(y), ( 8 ) 

dy 

where h{y) = dy' fo{y', y). 

The Thomas-Chao-Quigg relationship, D^iy) = 
ndNch/dy, can be solved explicitly in two extreme cases 
when either Mq = 0 or Mch = 0. When Mq = 0, we have 

/ y foo 

dy'giy') dy" g{y") = Kg{y) (9) 

-OO J y 

and the solution is given by 

giy) = ^sech^(|:) ( 10 ) 

which can be easily verihed using sech^a: = dtauhx/dx = 
1 — tanh^ a;. With k = 0(1), this form alone (basically 
the modihed Poschl-Teller potential) is much too sharp 
to describe a realistic dNch/dy. Eurthermore, it has no 
room for energy dependence once k is fixed. This is in 
contradiction with the dNch/dy spectrum in elementary 
collisions which shows no prominent central peak of a 
fixed width. Hence the Mq = 0 scenario can be excluded. 
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In the Mch = 0 limit, the Thomas-Chao-Quigg rela¬ 
tionship is 

/ y POO poo 

dy' dy''fo{y',y'') = K dy'fo{y',y) (11) 

-oo J y J — oo 

To solve for fo{y,y'), we make an ansatz 

/o(y,y') = ^(l/rei)F(r), (12) 


where ?/rei = y — y' and Y = [y + 2/0/2- The nor¬ 
malization conditions for R and F are dy R{y) = 
dy F{y) = 1. Equation E itTTll can then be solved by 
making change of variables to y^ei and Y. The solution 
is 


fo{y,y') = (■ 




(13) 


where the only restriction on F is that the integrals 
in Eq. cu converge and it reproduces the experimental 
dNch/dy. For details, see Appendix [Cl This form of 
correlation function is very reasonable as this is nothing 
but the distribution function of the decay products of a 
cluster whose rapidity is distributed according to F(Y). 

For small enough k, we should have 


F{y) 


1 dNch 
Nch dy 


(14) 


where Nch = 2 {Mq) is the total charge multiplicity. In 
this solution, it is clear that k is directly related to the 
correlation length A in the relative rapidity space of the 
pair y — y' as 


\ — 2k 


(15) 


Taking the correlation k ~ 1, it is easy to see that 
the charged particle spectrum will be a smooth distribu¬ 
tion with typical variation in rapidity of 2 k ~ 2. This is 
in good agreement with the pp collision results in FigH 
Further, the absence of a narrow peak in the central re¬ 
gion also excludes significant contribution from the un¬ 
correlated charged particles. 

If we have a finite observational window (—j/o, 2/o), the 
charge transfer fluctuations are given by: 


Du{y) 


(Mo) 

2 



dy"h{y',y") 


/ Vo pco 

dy' dy"U{y',y") 

■Vo 


/ y rvo 

dy' dy"h{y',y") 

-Vo '^y 


(16) 


where the bar in indicates that this quantity is mea¬ 
sured in a limited window. With a moderate yo, the 
Thomas-Chao-Quigg relationship Eq. o no longer holds 
even if the unrestricted charge transfer fluctuation sat¬ 
isfies it. However Du{y) is still sensitive to the charge 


correlation length. To have an idea how Du{y) behaves, 
we can use a flat F{Y) following Thomas and Quigg. If k 
is much smaller than the size of the total rapidity interval 
Ttnax and yo is not too close to Ycaa^l2, we get 


Du{y) 

{Nch) 


— I"- -h 

2yo [2 2 


2e cosh 



If we have a large //o/k, this becomes constant and we get 
back to the original Thomas-Quigg relationship Eq. 
With a finite j/o, this is a monotonically decreasing func¬ 
tion of y > 0 with the maximum at y = 0. 

The pseudo-rapidity distribution measured by RHIC 
experiments does show a plateau within — 2 < 77 < 2 and 
the rapidity distribution shows a similar plateau within 
— 1 < y < 1. Hence, if the correlation length (k) remains 
constant throughout the plateau, one would expect that 
Du{y) measured within the plateau should also have a 
maximum at y = 0. On the other hand, if one observes 
that Du{y) has a minimum at y = 0, then it is a signal 
that quite a different system (a QGP) is created near the 
central rapidity with a much smaller rapidity correlation 
length than the rest of the system. 

To examine the relationship between the charge trans¬ 
fer fluctuations and the charge correlation length, we 
must make sure that other effects such as the impact 
parameter fluctuations and the hadronic correlations do 
not mimic the effect we seek. To study the non-QGP ef¬ 
fects, we have analyzed 50, 000 HIJING events |^. The 
results are shown in FigEl Each point in the figure repre¬ 
sents 5% bin in the centrality measured by the number of 
charged particles within — 1 < y < 1. It is quite obvious 
from this figure that the charge transfer fluctuations do 
not vary with centrality. It is also obvious that 5„(y) in 
this case is a decreasing function of y. 

It is interesting to compare Eq. with the results 
from HIJING. For this, we used the top 15% central re¬ 
sults in FigISl and fitted them with En. dl 711 . The best 
fit gives K = 0.72 though the results from Eq. 03 are 
slightly flatter than the HIJING results. This discrep¬ 
ancy in shape is not unexpected because the HIJING 
dNch/dy is not well approximated by the flat F{Y). 
The full pseudo-rapidity space analysis of 50,000 HIJING 
minimum bias events are shown in Fig0 for 3 different 
centralities. The Thomas-Chao-Quigg relationship works 
quite well within this model for all centrality classes. In 
the rapidity region of y = 0 ~ 3, the value of k is in 
the range of 0.63 ~ 0.68. This is consistent with the 
experimental pp result. 

In the next two sections, we test our idea against two 
scenarios for heavy ion collisions. In the first scenario, 
the created system consists of a single species of neutral 
clusters which may be taken as hadronic clusters. In the 
second scenario, the created system contains a second 
component with a much smaller correlation length. 
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Nch'^Nmax 

FIG. 3: The ratio of the Charge transfer fluctuations to the 
number of charged particles Du{y )/is shown as a function 
of centrality A'ch/A^max- The observation window is fixed at 
(—1.0,1.0). These results are from a HIJING calculation for 
y/s = 200 GeV An-An collisions. Parameters for the HIJING 
calculations are dE/dx = 2 GeV/fm, pminijet = 2 GeV with 
the shadow and the quench flags on. The most central bin 
has the top 5 % of the events. All other bins have the size of 
10%. Total of 50,000 HIJING events are analyzed. 



FIG. 4: The results of analyzing simulated Au-Au events 
using HIJING. The lines are scaled pseudo-rapidity spectra 
for centrality classes 0 — 5 %, 25 — 35 % and 45 — 55% and the 
symbols are charge transfer flucuations for each class. 


III. SINGLE COMPONENT MODEL 

In this section, we consider the charge transfer fluctu¬ 
ations in a system which consists of only a single species 
of neutral clusters (presumably a hadronic matter). 

The Thomas-Chao-Quigg relationship was used to jus¬ 
tify the use of the neutral cluster model in pp collisions 
in the last section. For heavy ion collisions, as far as 
we know there has been no experimental investigation in 
this area. For this study, we take the HIJING simulation 
results shown in Fig 0 as an indication that the single 


component neutral cluster model is also a good hadronic 
model of AA collisions and explore the consequences. A 
case study with a few hadronic event generators is under 
way and will be reported elsewhere. 

With only single species of clusters, the joint proba¬ 
bility to have an event with particle rapidities ({j/a}) is 
simply 

Mo 

Piiya}) (18) 

fc=l 

All observables in this model depend only on the pair 
distribution function ^y~)- Hereafter, we will drop 

the subscript 0 from /o for brevity. This is pertinent to 
the heavy-ion collisions at RHIC where the net charge is 
almost zero in the central region. For completeness, we 
have also listed the formulae for unpaired net charges in 
the appendix. As discussed in the appendix, the contri¬ 
bution of the non-zero net charges to the charge transfer 
fluctuations is negligible as long as Q ^ Wh- 

Due to our assumption that the decay products are 
identical except for the electrical charges, the 2 -point 
function has the following symmetries 

= f{y~,y^), (i9a) 

f{y^,y~) = f{-y^,-y~)- (i9b) 

We will use these symmetry properties to simplify the 
equations for the charge transfer fluctuations and the net 
charge fluctuations. 

Guided by our discussion of the Thomas-Chao-Quigg 
relationship, we use a separable form of pair distribu¬ 
tion function f{y'^,y~) = F{Y)R{y^ei), where Y = 
{y'^ + y~)l‘^ and yrei = y~^ — y~ as in Eas. ltT^ and 
m- The function F{Y) is the rapidity distribution of 
the clusters and the function i?(?/rei) is the relative rapid¬ 
ity distribution of the produced pair. For the shape of 
F{Y), we use a Wood-Saxon form and a gaussian form 
here but more sophisticated forms are certainly possible. 
For i?(?/rei), the following two physically motivated forms 
were used in our study: 

R{r) = exp (-rV27') (20a) 

R{r) = exp(-|r|/ 7 ) ( 20 b) 

27 

The function in Ea. ll20bll is the same as Eq. m with 
7 = 2k. 

The combination of F{Y) and R{r) are not arbitrary. 
Once we fix the charge correlation length 7 for a partic¬ 
ular R{r), then the parameters for F(Y) are complete 
determined by the best fit to the experimental charged 
particle distributions |H 3 • To test the sensitivity to 
the different forms of the correlation function, we use 
the following 4 parameter sets in this section: The pa¬ 
rameter set labeled 1 uses a Wood-Saxon form of F{Y) 
and a gaussian i?(?/rei)- The set labeled 2 uses a Wood- 
Saxon F{Y) but i?(yrei) bas the exponential from. The 
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parameter set labeled 3 uses a gaussian F{Y) and also 
a gaussian Riy^ei)- The set labeled 4 uses a gaussian 
F{Y) but R{yre\) has the exponential form. These will 
be used for the net charge fluctuation analysis and the 
charge transfer fluctuation analysis. 

A. Net Charge Fluctuation 

STAR collaboration at RHIC has published their mea¬ 
surement of net charge fluctuations and concluded 
that their result is consistent with the hadronic gas ex¬ 
pectations. In this section, we use this data to fix the 
correlation length 7 . Since our previous fit to the HI- 
JING simulations gave 7 ^ 1.3 — 1.4, we will look for the 
correlation length within the range of 1 < 7 < 2. 

Within the observation window —yo < y < yo, the net 
charge and the charge transfer are given by 

u{y) = IQpiy) - Qsiy)] /2 
Q{yo) = [QF{y) + QB{y)] ( 21 ) 

where now Qriy) and Qsiy) are measured within the 
observational window. Notice that Q{yo) is actually in¬ 
dependent of y. In the limit y = yo, we have u = Q/2. 

In terms of the charge pair distribution function, the 
net charge fluctuation is given by 

SQ^yo) = {Q{yof) - {Q{yo)? 

/ -Vo rVo 

dy~ dy+f{y+,y~)(22) 

-00 —yo 

In deriving the above equation, we have made full use 
of the symmetry properties of the pair distribution func¬ 
tion, Eqs. IIH). Since the two particle distribution func¬ 
tion f{y~^, y~) is peaked at |y+ — y~\ =0, the net charge 
fluctuations in Ea. (l22ll measure the local correlations at 
around the edge of the observation window y ~ ±?/o- One 
can also see that as yo —>■ 00 , the net charge fluctuation 
vanishes as it must. 

The total number of charged particles in the given ra¬ 
pidity window (-2/o,2/o) is: 

/ Vo poo 

dy~ / dy+f{y''~,y~). (23) 

-yo J—oc 

The ratio of the net charge fluctuations to the total num¬ 
ber of charges in the rapidity window (—0.5, 0.5) was 
measured in the RHIC experiments, and the value for 
central collisions is around SQ^/Noh = 0.8 —0.9 after cor¬ 
recting for global charge conservation effect . The 

correction is made through a factor of 1/(1 — p) where p 
is the fraction of the charged particles included in the 
observation compared to the total number of charged 
particles produced. ,i,From this value, we can find the 
correlation length for charged pairs. The exact value de¬ 
pends slightly on the assumption on the shape of the pair 
correlation function R{r) in relative rapidity space and 
the charge center distribution F{Y). 



y 

FIG. 5: The charge conservation corrected ratio of the net 
charge fluctuations to the total number of charged particles 
5Q^(j/o)/Ach(j/o) is shown as a function of the charge corre¬ 
lation length 7 . We show the ratio at three different rapid¬ 
ity observation windows, \y\ < t/o, where yo = 0.5, 1.0 and 
1.5 for the three lines respectively. The parameters for the 
charge pair center distribution function F{y) are adjusted to 
fit the charged particle distribution data measured by PHO- 
BOS group at ^ = 130 GeV. 



Vo 

FIG. 6 : The corrected net charge fluctuation ratio 

(yo)/Nch{yo) is shown as a function of the observation 
window size t/o for different values of charge correlation 
lengths 7 = 1.25, 1.5, 1.75 and 2.0 respectively. The net 
charge fluctuations are decreasing functions of observation 
window. 


In FigEl we have plotted the ratio of the net 
charge fluctuations to the total number of charges 
(yo)/Nch{yo) as a function of the pair correlation 
length 7 . We only show here the charge fluctuation re¬ 
sults from the parameter set 2, where the charge center 
distribution F{Y) is Wood-Saxon form and the relative 
rapidity distribution between the pair R{r) is exponen¬ 
tial decay form. The different assumed forms of the pair 
distribution function / = F{y)R{r) have little effect on 
the net charge fluctuations and are not show here for 
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FIG. 7: The charge conservation corrected ratio of the net 
charge fluctuations to the total number of charged particles is 
plotted as a function of the charge correlation length 7 . Set 
1 refers to a Wood-Saxon form for the charge pair center dis¬ 
tribution F{y), and a gaussian form for the relative rapidity 
distribution between the pair R{r). Set 2 refers to a Wood- 
Saxon form for F{y) and an exponential decay form for R{r). 
Set 3 and 4 are for gaussian form of F{y) with gaussian or 
exponential decay form of R{r) respectively. This should be 
compared with data reported for ^/s = 130 GeV Au-Au colli¬ 
sions at RHIC. Here the observation window is (—0.5, 0.5) in 
rapidity. 


clarity. From this figure, one can conclude that the net 
charge fluctuations SQ^{yo) are strongly correlated with 
the charge correlation length 7 in this single component 
case. 

RHIC experiments can measure the net charge fluc¬ 
tuations as a function of the observational window size 
2 / 0 . We have plotted the corrected net charge fluctua¬ 
tions as a function yo in FigEl As can be seen in this 
figure, the net charge fluctuations always decrease when 
the observation window is enlarged. This is because the 
total number of charged particles included in the obser¬ 
vation window is increasing faster than the net charge 
fluctuations. The slope however is related to the charge 
correlation length. 

Using the corrected net charge fluctuation ratio data 
from STAR |^, we deduce that the charge correlation 
length is 

7Ril.5 (24) 

This deduced value is largely independent of the shape of 
the charge pair correlation function R{r) or the pair cen¬ 
ter distribution function F{Y) as shown in Fig[7| Notice 
that the inferred charge correlation length is consistent 
with the HIJING results where k = 7/2 « 0.7 in the pre¬ 
vious simple estimate within the rapidity window (— 1 , 1 ). 



FIG. 8 : The ratio of the charge transfer fluctuations Du{y) to 
the charged particle yield dN^h/dy in the observation window 
(—1.5,1.5) is plotted as a function of the charged particle 
correlation length 7 . The different lines represent the results 
for different forward-backward rapidity cut, y = 0.0, 0.3, 0.6, 
0.9 and 1.2 respectively. We only showed the results from 
parameter set 2 , while other parameter sets yield qualitatively 
similar results. This result is for the one component model. 


B. Charge Transfer Fluctuations 

We now apply our model to the charge transfer fluc¬ 
tuations using the same parameter sets for the pair dis¬ 
tribution function f{y'^, y~) as we have used in the last 
section. 

First, we need to express the the charge transfer 
fluctuations in terms of the pair distribution function 
y~ ). Making full use of the symmetries of the func¬ 
tion f{y'^,y~) given in Eg. 11911 . we can simplify Ea. 11611 
to 

Du{y) = ^-^^-{5QF{y)5QB[y)) 

= ^^^+2(Mo) r dy- rdy+f{y+,y-) 

^ -Vo 

(25) 

where we used Eq. ll^ and defined SX = X—{X). Writ¬ 
ten this way, it is clear that the charge transfer fluctua¬ 
tions depend on both the charge correlations at the edges 
of the observation window y ~ ij/o and at the forward- 
backward rapidity cut y. Erom this expression, one also 
sees that the lower limit of charge transfer fluctuations is 
{Du{y)}min = i5(5^(yo)/4. 

In EiglHl we have plotted the ratio of the charge 
transfer fluctuation Du{y) to the charged particle yield 
dNch/dy as a function of pair correlation length 7 . As the 
observed rapidity spectrum is nearly flat around midra¬ 
pidity, dividing by Nch or dNch/dy only affects the over¬ 
all scale. The value of the ratio strongly depends on 
the correlation length 7 . Also the ratio decreases as the 
forward-backward separation y increases at a fixed 7 . As 














mentioned before (c.f. Eq. C3), this decrease is due to the 
limited observation window and has nothing to do with 
the changing correlation length. This can be also shown 
in the following way. If dNch/dy does not vary signifi¬ 
cantly within the observational window, then F{Y) does 
not vary significantly within the observational window. 
In that case, 

1 dD(v) 

j— « -F(0) / dr R{r)<Q (26) 

(Mo) dy Jy^_y 

for y > 0. Therefore, if dNch/dy is flat and the system is 
composed of only one species of neutral clusters, Du{y) 
must be a decreasing function of y > 0. Conversely, if 
Du{y) is an increasing function of y > 0 while dNch/dy 
remains flat, it signals the existence of a new component. 
We turn to this possibility in the next section. 


IV. TWO COMPONENT MODEL 



y 

FIG. 9: The Charge yields are plotted as a function of ra¬ 
pidity in a two component system. The full line represent 
the total charged particle yield profile while the dashed and 
dotted lines represent the contributions from the hadron gas 
and the QGP phases in a typical two component model we 
used. The QGP size is ^ = 1.28 in this case. 


In a two component model, the full joint probability 
for the rapidities is given by 


All AA2 

p(.{yi}) = n Aivt^yj) n. y //)- (27) 

3 = 1 k=l 

where /i and /2 have different correlations lengths 71 > 
72 - The pair distribution functions are again taken as 
the separable form: fi{y~^,y~) = Fi{Y)Ri{r) where Y = 
( 2 /^ +2/~)/2 and r = y^ —y~ . Since the fluctuations add 
in quadrature, the net charge fluctuations and the charge 
transfer fluctuations are just the sum of contributions 
from the two components, fi and /2 respectively. 

Physically, the two component model is motivated by 
the fact that the hot and dense matter produced in the 
relativistic heavy-ion reaction is not necessarily homoge¬ 
neous in rapidity space as explained before. If the de- 
confined QGP phase did exist during the early stage of 
heavy-ion reaction, it is highly possible that the QGP 
phase coexisted with the hadron gas phase. A simple 
situation would be a phase separation between the QGP 
phase and the hadron gas phase. This could produce 
signals that are specific to a phase coexistence scenario. 
In the case of the charge transfer fluctuations, a QGP 
and hadron gas phase separation could be measured and 
mapped into the charge correlations in the relative rapid¬ 
ity space between the pair of particles. In the following, 
we will refer to the short correlation part as a ‘QGP’ and 
the long correlation part as a ‘hadron gas’ (HG). 

In our simple two component model, we assume the 
two components have different correlation lengths and 
the R{r) functions are either taken as an exponential 
form or a gaussian form, with the corresponding corre¬ 
lation lengths satisfying 71 > 72 . Here 71 = 7 hg is the 
rapidity correlation length of the hadronic part (labeled 
‘1’) and 72 = 7 QGP is that of the QGP part (labeled 


‘2’). We let the cluster distribution functions for the two 
components be: 


(Mi)Fi(y) 


(M 2 )F 2 (y) 


_Cl_ 

I -b exp [(|y| - cro)/ao] 


C292iy) ■ 


C29i{y), 


(28a) 

(28b) 


Where ci is a normalization factor and C 2 is the strength 
of the QGP phase. To be physically consistent, the value 
of C 2 is adjusted so that the function Fi is always pos¬ 
itive. We then demand that f Rigi = f i ?252 so that 
dNch/dy is independent of the choices of gi and 92- The 
parameters (Tq and oq are chosen to fit the experimental 
data. 

For a gaussian Ri{r), it is convenient to take the func¬ 
tional forms of gi{y) and 92 ( 9 ) to be also gaussian. To 
satisfy / R 191 = f R 292 , the widths of these gaussians 
should be related as follows = 7 ^/4 -b = 7 !/4 + cr| 
where di is the width of gi{y) and ^ is the width of the 
QGP part of dNch/dy. For an exponential Ri{r), deter¬ 
mining the forms of gi and 92 is a little more complicated 
than the gaussian case. However, since this form of the 
correlation function satisfies the Thomas-Chao-Quigg re¬ 
lationship, we will use mainly the exponential form here¬ 
after. As in the one component model, both the net 
charge fluctuation and the charge transfer fluctuation re¬ 
sults are not very sensitive to the particular choice of the 
functional form the two particle distribution function. 

With the exponential Ri{r), the charge center distribu¬ 
tion functions gi{y) and 92 ( 9 ) can not be both gaussian 
anymore. It is convenient to select 9 i{y) (the hadronic 
part) to have a gaussian form. Using the fact that the 
function Ri{r) is in fact a Green function of the differ¬ 
ential operator {d/dr)^ — l/yf, the function 92 {y) (the 
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QGP part) can be obtained as 



9i{y) + 


1 - 



(29) 


where p{y) = J Rigi = f i? 2 g 2 is 


p(y) 



/ x + y \ _ 

V 2 y) 271 


7i J 


(30) 


The function p{y) is proportional to dNch/dy of the QGP 
part. The width ^ of the QGP part is determined by the 
width of p{y). 

In all, we have 8 parameters here. We have ci, tro, oo 
and C 2 explicitly appearing in Eq. (P|). The parameter 
(Ji is the width of the function gi. The parameter ^ is 
the width of the function p(y) and is connected with the 
parameter a-i through Eq. CT . We also have two charge 
correlation lengths 71 and 72 with the condition 71 > 72 . 

Among the 8 parameters, 5 are fixed in the follow¬ 
ing way. Since we are only interested in the ratio 
Du{y)/{dN/dy), the value of parameter ci is irrelevant 
and we just fix it to be 1. The parameters oq and ao are 
fixed by requiring that the resulting dN/dy shape de¬ 
scribes results from RHIC experiments. The parameter 
C 2 is always chosen to be the maximum possible value for 
the condition Fi{y) > 0 given all other parameters. The 
parameter cti is determined by the parameter 

The three parameters we are going to vary in the fol¬ 
lowing are then 71 , 72 and An example of total charged 
particle spectrum with the respective hadron gas and 
QGP contributions are plotted in Fig El There is a sub¬ 
stantial presence of QGP around midrapidity, but it be¬ 
comes less prominent for |y| > 


A. Net Charge Fluctuation 

As in the one component model, we first consider the 
net charge fluctuations to further fix our parameters. 
When there are two distinct species of neutral clusters, 
the net charge fluctuation within the rapidity interval 
(-yo, 2 /o) is given by 

/ -Vo fVo 

dy~ / dy+fi{y+,y~) 

-00 J-yo 

/ -Vo rVo 

dy~ / dy'^f 2 {y'^,y~) (31) 

-00 —Vo 

and the total charged multiplicity is given by 

/ Vo roo 

dy~ / dy~'~ fi{y~,y+) 

-yo J—00 

/ yo roc 

dy~ / dy~'~ f 2 {y~,y~'~) (32) 

-yo J — OQ 

The ratios 5Q'^{jjo)/Nc\i{yo) corrected as before by a 
factor of 1/(1 — p) are plotted as a function of the ob¬ 
servation box size j/o in FiglTUI The lowest solid line 



FIG. 10: The ratios of the corrected net Charge fluctuations 
to the total number of charges 5Q^[yo)/Nch{yo) in the rapid¬ 
ity observation window {—yo,yo) are plotted as a function of 
rapidity yo for the one and two component models. The line 
for the one component model is labelled by the charge corre¬ 
lation length 7 = 1.44, and the lines for the two component 
models are labelled by the pair of charge correlation lengths 
in the two components (7hg,7qgp/7hg). The corrected net 
charge fluctuations are fixed 5Q^(0.5)/Ach(0.5) = 0.85 for all 
lines. 


is for the one component model with charge correla¬ 
tion length 7 = 1.55 as obtained in the last section. 
For the two component model, 4 different choices with 
7 i = 7 hg = 1.75,2.0 and 72/71 = 7 QGp/ 7 hg = 0.25,0.5 
are shown. Since the corrected net charge fluctuations 
at 2/0 = 0.5 is about 0.8 — 0.9 [3 |3^, the QGP width 
parameter ^ are chosen in such a way that all net charge 
fluctuations have SQ'^(yo)/Nch{yo) = 0.85 at yo = 0.5. 
For instanec, for the parameter set 7 hg = 1-75 and 
7 QGP = 7 hg /4 = 0.44, the width of the QGP part ^ 
turned out to be ^ = 1.28. 

The net charge fluctuations as a function of rapidity 
is flatter in the two component model than in the one 
component model. This is because the two component 
results interpolate between the one component results 
with 7 = 7QGP and 7 = 7hg- Since the behaviors of 
the net charge fluctuations for the one component model 
and the two component model are clearly distinct, one is 
tempted to argue that the flat (5Q^(?/o)/Vh(2/o) itself is 
an indication of a second phase. (A similar idea was sug¬ 
gested in Ref. [ 5^.1 Unfortunately, totally uncorrelated 
system also has a flat SQ^(yo)/Nch{yo) when corrected 
for the effect of total charge conservation. 

In addition, the net charge fluctuation 
{yo)/Nch{yo) is constrained by the fact that in 
the limit yo 0, we should get the Poisson limit 

{yo)/Nch{yo) —> 1- This puts a constrain on the 
sensitivity of net charge fluctuations to the QGP phase. 
In our two component model, the QGP phase is located 
mostly around midrapidity, and the presence of QGP is 
reduced at larger rapidities. Hence in order to observe 
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the QGP, we need to have yo ~ 0. But because of the 
limiting value at = 0 , the net charge fluctuations 
actually have a reduced sensitivity to the QGP phase. 
The charge transfer fluctuation, which we now turn our 
attention to, does not have these limitations. 


B. Charge Transfer Fluctuation 


The charge transfer fluctuation Du{y) is qualitatively 
different than the net charge fluctuation SQ^{yo)- As 
will be shown in this section, the charge transfer fluctu¬ 
ations is capable of distinguishing the two phases of our 
two component model and hence can be used as a signal 
for the QGP phase. In our model, the charge transfer 
fluctuation Du{y) is given by 

Du{y) = + f dy- f dy+fi{y+,y-) 

/ y rVo 

dy~ / d?/+/ 2 (y+,y“) (33) 

-Vo Jy 

while the final particle spectrum is 

^ = 2(Mi)/ii(j/) + 2(M2)h2(y). (34) 

ay 

where hi{y) = J^^dx fi{x,y). 

For small 7 i, it is easy to show that 

Du{y) ~ const-f 7 i(Ml)Fi(y) -f 72 (M 2 )F 2 (y) (35) 


while the rapidity spectrum becomes 

^ ^ (Mi)Fi(j/) + (M2)F2 (j/) . (36) 

dy 

Hence, changes in the charge transfer fluctuation com¬ 
pared to the charged particle spectrum reflect changes 
in the concentration of the two components and/or the 
change in the mean correlation length. When the corre¬ 
lation lengths are not very small, Du{y) is given in terms 
of the convolution of Fi{y) and ^ 2 ( 2 /) with the corre¬ 
sponding relative distribution Ri{r) and R 2 {r). Unless 
7 ’s are very large, the ratio Du{y)/{dNch/dy) should still 
be sensitive to the changes in the composition. 

If the net charge fluctuation SQ^{yo) is sizable, its pres¬ 
ence can reduce the sensitivity of the ratio 


K{y) = 


Dujy) 

{dN^Jdy) 


(37) 


to the changing composition since dNch/dy abruptly de¬ 
creases beyond the central plateau. However, since this 
^Q^iVo) is the uncorrected net charge fluctuation, it is 
easy to measure and subtract it from Du- In this case, 
the relevant ratio becomes 


K{y) 


Dujy) - SQ‘^(yo)/4: 
(diVeh/dy) 


(38) 
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FIG. 11: The ratio of the Charge transfer fluctuations to 
the total number of charges Du{y)/{dNch/dy) in the rapid¬ 
ity observation window (—1.5, 1.5) is plotted as a function 
of the forward-backward separation cut y in the one com¬ 
ponent and two component models. For the one component 
model, the charge correlation length is fixed at 7 = 1.55. For 
the two component models, the pair of correlation lengths 
(7hg, 7qgp/7hg) are labelled on each line. For the mid¬ 
dle two lines, the QGP sizes ^ are fixed by requiring the 
net charge fluctuations to be constant (5Q^(0.5)/A^ch(0.5) = 
0.85. The two lines with the same charge correlation lengths 
(7hg,7qgp/7hg) = (1.75,0.25) have different QGP size. 


If dNch/dy is flat within the observational window, this is 
of course not necessary as SQ^{yo) term just adds a con¬ 
stant. Also in the large yo limit, SQ'^ijjo) 0 due to the 
overall charge conservation and hence this modification 
is not necessary. 

The quantity K(y) within —1.5 < y < 1.5 are plotted 
in Figtm as a function of the forward-backward rapidity 
separation y in the one and two component models. As 
our dNch/dy is almost flat within this window, we don’t 
need to subtract SQ‘^{yo) part. The shape for the one 
component model is completely fixed by the charge cor¬ 
relation length 7 = 1.55 as before, and it is a decreasing 
function of y. For the dashed and the dotted lines, we 
use the same parameters as obtained in the last section 
based on the experimentally observed net charge fluctu¬ 
ations. Even though the the one and two component 
cases have a common net charge fluctuation at yo = 0.5, 
the charge transfer fluctuation patterns are quite differ¬ 
ent: The most prominent feature for the two component 
model is the appearance of the minimum for k(y) at y = 0 
for 7 QGP < while the single component case al¬ 

ways has a maximum at y = 0 . 

The minimum appears at midrapidity because that 
is where the QGP component is concentrated. As y 
increases, the fraction of the QGP matter decreases. 
Hence, K(y) increases as a function of y > 0. The point 
where the slope of ^(y) changes sign must is directly re¬ 
lated to the width of the QGP component. Unfortu¬ 
nately, when the size of the observation window and the 
width of the QGP component are similar, the sensitivity 
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FIG. 12: The charge transfer fluctuations Du{y)/(dNch/dy) 
are shown as a function of the forward-backward rapidity cut 
y in the two component system. The QGP size is indicated 
by the parameter We have also plotted the corresponding 
hadron gas charge profiles, dNch/dy{HG), for the two sets of 
calculations. 


to the size of QGP is partially lost because R,{y) must 
decrease as y approaches the edge. To measure the size 
of the QGP component well, one needs to have j/o 
This prompts us to extend the rapidity window of our 
observation in the two component calculations. 

In Fig ll2l we plot K{y) with two different ^’s as a func¬ 
tion of y. The observation window used is (—6.0,6.0) 
which are large compared to the size of the QGP compo¬ 
nent. One can conclude that the width of the depression 
around midrapidity does reflect the width of the QGP 
component. Nearly flat K{y) for y > |^| is expected since 
in this region we should recover the single component 
result. The decrease near the edge is again due to the 
finite window size. Also shown are the shapes of the 
HG contribution to the rapidity spectrum as a reference. 
The charge transfer fluctuations with the same two ^’s as 
above are also shown in Figim In the case of a more lim¬ 
ited window, the sensitivity to the QGP size is reduced 
due to the fact that the size of yo is in fact about the 
same as the width of the QGP part 

The charge transfer fluctuation Du{y) is completely 
different from the net charge fluctuation SQ‘^{yo) as far 
as the observation window size effect is concerned. As 
discussed before, having a larger observation window can¬ 
not increase the sensitivity of the net charge fluctuations 
to the QGP phase when it is confined to a small region 
around midrapidity. However, for charge transfer fluctu¬ 
ations, having a large observation window increases the 
sensitivity since the window now encompasses more of 
the QGP part. Enlarging the observation window also re¬ 
duces the edge effect increasing the sensitivity even more. 
An additional advantage for the charge transfer fluctu¬ 
ations is that there is no the global charge conservation 
correction unlike the net charge transfer fluctuations. 

For a further reference, we show the result of analyz¬ 



FIG. 13: A plot of ratios {Du{r}) — SQ'^/4:)/{dNch/dr]) for 
3 different centrality classes using 50,000 minimum bias HI- 
JING events for An-An collisions at = 200 GeV. 


ing 50,000 minimum bias HIJING events for Au-Au col¬ 
lisions at y/s = 200 GeV in Fig^] for 3 different cen¬ 
trality classes. Again the analysis is carried out in the 
pseudo-rapidity space. The data points used are the same 
as in Fig0 except that the charge transfer fluctuations 
are corrected for the overall net charge fluctuations as in 
Ea. (l38ll . The fact that SQ^ ^ 0 even in the full phase 
space is due to the spectators. Since the net charge car¬ 
ried by spectator nucleons fluctuates, so does the net 
charge of the produced particles. In the usual way of 
characterizing the centrality classes (by Vch or Et), this 
is unavoidable. It should be observed that for all central¬ 
ity classes, the ratio K{ri) = {Du{ri) — SQ^/A)/{dNch/dy) 
is essentially flat. It is somewhat surprising the Thomas- 
Ghao-Quigg relationship works well for HIJING events 
given the fact that no ‘cluster’ appears explicitly within 
HIJING. 


C. Forward-Backward Multiplicity Correlation 

In all the above considerations, the key points are: (i) 
there are primordial clusters that produce multiple parti¬ 
cles, (ii) a local cut separates the phase space into two re¬ 
gions, (iii) one can define observables that only count the 
primordial clusters that have their decay products sepa¬ 
rated by the local cut (c.f. FiglU). These points imply 
that such observables are sensitive to the local proper¬ 
ties around the cut. Hence, if the nature of the ‘clusters’ 
changes in different regions of the phase space, then these 
observables can detect the changes. 

In some experiments, such as PHOBOS at RHIC, 
charge states of the produced particles cannot be deter¬ 
mined. In this case, charge transfer fluctuation cannot be 
used. However, since the essence of the current method 
is to have a cut that separates produced particles, just 
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measuring the forward-backward multiplicity correlation 
w{f]) = {Nf{v)Nb{v)) - {Nf{v)){Nb{v)), (39) 


may be enough to detect the change in the correlation 
length. We explore this idea in the following. Here NF{ri) 
is the charged multiplicity in the region forward of the 
pseudo-rapidity ij and NBirf) is the charged multiplicity 
in the region backward of rj. 

In this section, we switch to the pseudo-rapidity ry since 
without particle identification one cannot determine the 
rapidity y. However, in the current formulation of the 
problem, the only change this switch introduces is that 
instead of rapidity correlation function fi{y,y') we have 
the pseudo-rapidity correlation function fi(ri,ri'). 

Using the two component model Eq. 123, it is easy to 
show that 

w{r]) = 4:[{SM^) - (Ml)) / dr]hi{r]) / dy hi{r]) 

Jr) J-r}o 

fVo rf] 

-h4 ((JMa) - (M 2 )) / dyh2{y) / dyh2{y) 
Jr) J—rjo 

fVo fV 

+2{Mi) / dy' / dy” fi{y',y") 

Jr) J-Vo 

fVo rv 

+2{M2) dy' dy" f 2 {ri',y") (40) 

Jr} J-r]o 

where the correlation functions fi and /2 are functions 
of pseudo-rapidities and we defined 

/ OO 

dv'Mv'^v) (41) 

-00 


The non-trivial part of this expression is essentially the 
same as the charge transfer fluctuations. The sensitivity 
of this observable to the changes in the pseudo-rapidity 
correlation length depends crucially on the size of the 
first two terms containing (SMf) — {Mi). If the number 
fluctuations of the clusters obey Poisson statistics, then 
these two terms vanish. In that case, w{y) is as sen¬ 
sitive as the charge transfer fluctuation to the presence 
of the second phase. In the limit where there is only a 
single species of clusters and also ryo ^ 00 , we have an 
additional Thomas-Chao-Quigg relationship 


w{y) 


K 


dNch 

dy 


(42) 


with a constant k, provided that {5M‘^) = (M). 

The ratio of the forward-backward multiplicity correla¬ 
tion w{y) to the charged particle yield dNch/dy is plotted 
in Fig- ITU When the clusters are distributed according to 
Poisson distributions, the Thomas-Chao-Quigg relation¬ 
ship holds for a single component model and the ratio of 
w{y)/(dNch/dy) is flat in the central region (The solid 
line in FigtTHl. For a two component model with Pois¬ 
son statistics (long dashed line), we see a minimum at 
midrapidity just as in the charge transfer fluctuation. 



y 

FIG. 14: The ratio of the forward-backward multiplicity cor¬ 
relation w{ri) to the charged particle yield dNch/dy is plot¬ 
ted is plotted as a function of y for one component and two 
component models. The fluctuations of the total number 
of charged particles are also indicated for the corresponding 
lines. For Poisson distribution of Mi (i = 1 or 2), the devia¬ 
tion from Poisson is 0%, that is, (SMf) — (Mi) = 0. We have 
also included the results for 15% deviations. 


When the statistics deviates from Poisson, the factor 
(SM/) — (Mi) is non-zero. Then the first two terms in 
Eq.® contributes. These terms decrease with y and 
hence partially compensates the rising part due to the 
correlation change. However, the qualitative trend of the 
forward-backward multiplicity correlation still remains 
valid: The increasing segment is still present in the two 
component model and is an indication of the presence of 
QGP. 

For heavy particles originating from a thermally equili¬ 
brated system, the multiplicity fluctuation should follow 
the Poisson statistics. For light particles, (SM'^) can de¬ 
viate up to 15 % from the Poisson value. As an estimate, 
in Fig. d we show our results with (SM/) — (Mi) = 
0.15 (Mi). One can still see a clear dip near midrapidity. 

Coincidentally, the PHOBOS group has also measured 
a variation of the signal we proposed here (see 
The difference between the forward-backward multiplic¬ 
ity correlation w{y) and the “charged particle multiplic¬ 
ity fluctuations” cr{C) used by the PHOBOS group is 
subtle but results in quite different sensitivity. The signal 
cr(C) in measures the correlations between rapidity 
regions of {—y — Ay/2, —y + Ay/2) and (y — Ay/2, y -b 
Ay/2), with each region covering the same rapidity win¬ 
dow of Ary = 0.5, 1.0 and 1.5. Typically these two ra¬ 
pidity regions do not have a common edge. In our case, 
the two pseudo-rapidity regions are {—you'll) iv^Vo)- 
They share a common edge at y, but the two regions 
are generally of different size. When y = Ay/2, a(C) 
is the same as the forward-backward multiplicity corre¬ 
lation r(;(0). However, as we have shown in this paper, 
the single point in the fluctuation measurement can not 
distinguish between one component and two component 
models as the charge correlation length can be adjusted 
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to fit this single point. One must measure w{r]) as a 
function of rj to get full information on possible phase 
change. 

V. SUMMARY 

In this paper, we proposed the charge transfer fluc¬ 
tuations as a signal for the QGP phase of matter. The 
essence of our argument is very simple. Suppose there are 
strong unlike-sign correlations in the underlying system, 
then the charges are locally conserved. With a separating 
wall in the local region where we want to explore, the cor¬ 
relations (or fluctuations) of charges across this wall are 
sensitive to the local charge correlation length. The pairs 
created far away from the separating wall contribute lit¬ 
tle to the fluctuations since they have little chance to be 
separated by the wall. This is essentially the idea behind 
the charge transfer fluctuations proposed for earlier pp 
collisions. 

Since the charge transfer fluctuation Du{y) is a lo¬ 
cal measure of the unlike-sign correlation length, it can 
be used to detect the presence of a second phase in 
AA collisions. Quite generally, one can say that the 
presence of a local minimum at y = 0 for the ratio 
K{y) = Du{y)/{dNch/dy) is a signal for the second (pre¬ 
sumably QGP) phase. The appearance of this minimum 
is due to the facts that (i) a QGP phase should ap¬ 
pear around midrapidity where the density is the highest 
(ii) hadrons coming out of a QGP phase should have a 
markedly short unlike-sign correlation length compared 
to that of the hadronic matter. The size of the depression 
near midrapidity in turn contains information on the size 
of the QGP phase. If one has a large observation win¬ 
dow, the extent of the second phase in the rapidity space 
can be in fact estimated by the width of the dip. 

Extending the idea of charge transfer fluctuations, we 
also proposed the forward-backward multiplicity correla¬ 
tion as a possible signal for the presence of a QGP. A 
case study with an embedded QGP component in a few 
hadronic event generators is under way. 

APPENDIX A: EFFECT OF NET POSITIVE 
CHARGES 

For the uncorrelated charged particles in the system, 
keeping the total charge constant in the system, we find 
the charge transfer fluctuations to be 

Du{y) = ^ r dy'g{y') 

^ J-Vo 

(Al) 

In deriving the above result. We have assumed that the 
positive and negative charges have the same normalized 


distribution, g{y). Otherwise, we have to count the con¬ 
tributions from both positive and negative charges sepa¬ 
rately in Eq. and additionally there will be an extra 
term corresponding to the difference of the positive and 
negative charge forward-backward asymmetry. 

The first term in Eq. m comes from the overall fluc¬ 
tuations of charged particle number in the observation 
window and its value is the same as in the net charge 
fluctuations case (except a factor of 4). The second term 
is due to the forward-backward asymmetry nature in the 
charge transfer definition. In the limit that the obser¬ 
vation window is sufficient large, Eq. mi) reduces to the 
second integral in Eq. Q. 

In the case that the system has both uncorrelated 
charges and correlated charges, we only need to add the 
result in Eq. EH to the previous result for correlated 
charges, Eq. ll^ . This will not change the quantitative 
features of the charge transfer fluctuations as a function 
of rapidity. The first term in Eq. EH is independent of 
y and will not affect any of our discussion except adding 
a constant. The second term is a decreasing function of 
y. So, the contribution from uncorrelated charges is still 
a decreasing function of forward-backward rapidity cut 
y. This is in line with the results for correlated charges. 
The decreasing trend of the charge transfer fluctuations 
as a function of the forward-backward rapidity cut y in 
uniform system is unchanged by this additional contribu¬ 
tion. The existence of an increasing segment will still be 
a signal of a second phase with smaller charge correlation 
length. 

Since in realistic heavy-ion collisions at RHIC ener¬ 
gies, most of positive and negative charges are created 
together with an opposite charge to conserve the net total 
charges, we can safely assume that uncorrelated charges 
are rare. The net positively charged particles originating 
from the projectile and the targets is only a small frac¬ 
tion of the total number of charged particles in RHIC 
energy heavy-ion collisions. For this reason, the correc¬ 
tions from uncorrelated charges are ignored in the most of 
this study. The qualitative features of the charge transfer 
fluctuations will not be sensitive to this correction term. 

We can make a simple estimate of the corrections from 
these uncorrelated charges assuming the uncorrelated 
charges are from the protons in the initial collision sys¬ 
tem. The maximum of the second term in Eq. EH scales 
as p^M+ /Mq, where p is the fraction of observed uncorre¬ 
lated charges to the total number of uncorrelated charges. 
In RHIC energy heavy-ion reactions, M+ /Mq ~ 0.04 and 
p is typically around 5% in the central region. Indeed, 
the corrections from uncorrelated charges is quite small, 
of order 10“"^. The correction to the net charge fluc¬ 
tuations from the uncorrelated charges are on the same 
order of magnitude as the corrections to the charge trans¬ 
fer fluctuations. The net charge fluctuations Dc{yo) and 
the total number of charges A^ch(yo) both acquires addi¬ 
tional terms and they are both equal to 4 times the first 
term in Eq. inj. 
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APPENDIX B: NONZERO CHARGE TRANSFER 
CASE 

When the charge transfer u{y) in Eq.lj^J does not aver¬ 
age to zero, the charge transfer fluctuations will acquire 
additional terms that are quadratic to the average charge 
transfer. 

In the neutral cluster model, the full result for the 
charge transfer fluctuations is: 

Duiy) = ^{WL + WR + 2Wy) 

+ ((<5M2) - (Mo)) (Bl) 

The weights for left and right edges observation window 
(—j/o, Vo) and for the forward-backward rapidity cut y are 
defined as: 

/ -Vo rVo 

dy' dy”Uy’,y”) 

-cxD J —yo 

/ Vo roQ 

dy' / dy" fo{y',y'') 

-Vo Jyo 

/ y rVo 

dy' dy"foiy\y"). (B2) 

-Vo 'Jy 

The last term in Eq. m stems from the nonzero average 
charge transfer in the system. An estimate would give 
{u{y)) < 10 and {SMq) — (Mq) ^ O.IMq, and the error 
from neglecting this nonzero average charge transfer is 
typically less than 10“^. 


We can now Taylor-expand both dYF{Y) and 

F{z + r/2) with respect to r and get the following re¬ 
lationship between the moments of g{r) 

F 2 n +1 /R 2 n = 2 K {2n + 1) (C4) 

where 

poo 

Rs= drg(r)F (C5) 

Jo 

and we used the fact that g{r) is an even function. 

Note that 

poo 

/ =2”+iK”+in! (C6) 

^0 

so that 

^ 2^^+'^K^^+‘^{2n+iy. 

Jp°° dx 2^"+i«;^"+i(2n)! 

= 2K{2n+l) (C7) 

Therefore 

g{r) = C exp(-|r|/2«:) (C8) 

where C is a normalization constant. Hence, the solution 
of Eq. (tnil is given by 


Mx,y) = Af exp{-\x - y\/2K) F{{x + y)/2) (C9) 

APPENDIX C: SOLUTION OF EO.dTTll 

where F{Y) can be quite arbitrary as long as its deriva- 
The Thomas-Chao-Quigg equation for neutral cluster tives are all finite and the integrals in Eq. are well 
distribution function is given by defined. 



dxfo{x,y) 


= K dyfo{y,z) 


(Cl) 


We make an ansatz: 
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Mx,y)=g{r)F{Y) (C2) 
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poc pz+rj2 poo 

/ dr g{r) / dY F(Y) = k / dr g{r) F(z + r/2) 
Jo Jz — rj2 J— oo 
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